"Linearized" Dynamical Mean-Field Theory for the Mott-Hubbard transition 
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The Mott-Hubbard metal-insulator transition is studied within a simplified version of the Dynamical 
Mean-Field Theory (DMFT) in which the coupling between the impurity level and the conduction 
band is approximated by a single pole at the Fermi energy. In this approach, the DMFT equations 
are linearized, and the value for the critical Coulomb repulsion U c can be calculated analytically. For 
the symmetric single-band Hubbard model at zero temperature, the critical value is found to be given 
by 6 times the square root of the second moment of the free (U = 0) density of states. This result 
is in good agreement with the numerical value obtained from the Projective Selfconsistent Method 
and recent Numerical Renormalization Group calculations for the Bethe and the hypercubic lattice 
in infinite dimensions. The generalization to more complicated lattices is discussed. The "linearized 
DMFT" yields plausible results for the complete geometry dependence of the critical interaction. 

PACS numbers: 71.10Fd, 71.27. +a, 71.30.+h 



I. INTRODUCTION 

The correlation-induced transition from a paramag- 
netic metal to a paramagnetic insulator (the Mott- 
Hubbard transition [lj§]) has been intensively studied 
within the single-band Hubbard model [0-p|: 
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The model describes conduction electrons with spin a on 
a lattice with nearest-neighbor hopping matrix element 
tij and a local Coulomb repulsion U. One of the first 
approaches to describe the metal-insulator transition in 
the half-filled Hubbard model has been the Hubbard-Hl 
approximation |(|. The alloy-analogy solution predicts 
a splitting of the density of states in upper and lower 
Hubbard bands for large values of U. On decreasing U, 
the insulator-to-metal transition occurs when the Hub- 
bard bands start to overlap. The critical interaction is 
approximately given by the free bandwidth: U c ss W. 
The Hubbard-III approximation, however, fails to de- 
scribe the Fermi-liquid properties in the metallic phase. 

Later, the Mott-Hubbard transition has been described 
within the Gutzwiller variational approach by Brinkman 
and Rice |Q. Starting from the metallic side, the transi- 
tion is marked by a diverging effective mass. The critical 
interaction is found to be U c — — 16eo where eo is the ki- 
netic energy of the half- filled band per particle for U = 0. 
The Brinkman-Rice approach, however, fails to describe 
the insulating phase above U c . 

A Dynamical Mean Field Theory (DMFT), which be- 
comes exact in the limit of infinite spatial dimensions, 
has been developed for the Hubbard model fls|-[To||. The 
DMFT is able to yield a consistent description of the 
metallic Fermi liquid for weak coupling as well as of the 
Mott-Hubbard insulator for strong coupling. In practice, 
however, the solution of the mean-field equations is by no 



means a trivial task. In particular, for U i— > U c problems 
may arise since the "mean field" A(u>) which has to be 
determined self-consistently, develops a strong frequency 
dependence on a vanishingly small energy scale. 

The first calculations for temperature T = were per- 
formed using the so-called Iterated Perturbation Theory 
(IPT) @. Within the IPT the highly correlated Fermi 
liquid for U U c is characterized by a narrow quasipar- 
ticle peak that is well isolated from the Hubbard bands. 
As a consequence the insulating gap appears to open 
discontinuously at the metal-insulator transition. These 
characteristics of the transition have been questioned by 
various authors pd|-p^| so that the issue of the metal- 
insulator transition for T — (and also for T > jl(|[l4) ) 
cannot be regarded as settled at the moment. 

Qualitatively, the IPT scenario for T = is cor- 
roborated by recent non-pcrturbativc calculations using 
the Numerical Renormalization Group (NRG) method 
pl|,p6| . However, the critical value for the transition is 
found to be significantly lower as compared to the IPT 
result. On the other hand, the NRG value for U c is in 
remarkable agreement with the result of the Projective 
Self-consistent Method (PSCM) @. 

The value of the critical interaction for the Mott- 
Hubbard transition is of great interest. Using the meth- 
ods mentioned above, an approximate calculation of U c 
is possible. This, however, represents a rather difficult 
numerical problem, the solution of which still depends 
on the approximation used. Within the framework of 
DMFT, an exact analytical result for the precise value 
of U c is still missing. Even an approximate analytical 
expression is not available up to now. 

With the present paper we propose a simplified treat- 
ment of the mean-field equations ("linearized" DMFT) 
which allows to obtain an explicit expression for U c at 
zero temperature. A fully numerical treatment of the 
DMFT would leave us with a mere number for U c and 
would hardly show up the characteristic trends for dif- 
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ferent geometries unless a large number of cases were 
studied. Contrary, the linearized DMFT is able to yield 
at once the complete geometry dependence of U c . In our 
opinion this outweights the necessity for further approx- 
imations. 

The main idea of the linearized DMFT is to approxi- 
mate the hybridization function for the coupling between 
the impurity level and the conduction band by a single 
pole. This is detailed in Sec. II. The reliability of the new 
approach is estimated by comparing the analytical results 
for the Bethe and the hypercubic lattice in infinite dimen- 
sions with the available numerical values from the PSCM 
and the NRG in Sec. Ill A. A satisfactory agreement is 
found. In the following we then demonstrate the predic- 
tive power of the approach. The geometry dependence 
of U c is derived for a number of more complicated lattice 
structures: inhomogeneous Bethe lattices (III B and C) 
and hypercubic films in infinite and finite dimensions (III 
D). Finally, in Sec. Ill E we discuss a first correction be- 
yond the linearized theory. Sec. IV summarizes the main 
results. 

II. LINEARIZED DYNAMICAL MEAN-FIELD 
THEORY 

A characteristic feature of the metal-insulator transi- 
tion is that the quasiparticle peak appears to be isolated 
from the upper and the lower Hubbard band for U > U c 
and T = 0. Whether or not there is a real gap, i. c. 
zero spectral weight between the quasiparticle peak and 
the Hubbard bands, is difficult to decide with any nu- 
merical method but not very important for the present 
approach. Essentially, our approach is based on two ap- 
proximations: 

(i) We assume that in the limit U i— ► U c the influ- 
ence of the high-energy Hubbard bands on the low-energy 
(quasiparticle) peak is negligible. This can be specified 
as follows: Within the DMFT the Hubbard model is 
sclf-consistently mapped onto a single-impurity Ander- 
son model (SIAM). In the effective SIAM, we divide the 
conduction-electron degrees of freedom in a high-energy 
part -f/high (the Hubbard bands) and a low-energy part 
fliow (the quasiparticle peak). The Hamiltonian of the 
effective SIAM is then written as 

#SIAM = -ffhigh + -^high-imp + -^imp + #low + ^low-imp , 

(2) 

where the coupling of the impurity -ffj mp = e dd\d a + 
Ud^d-\d\di to the high- (low-)energy part is denoted as 
H high-imp (^low-imp)- The first approximation is then to 
neglect the terms if high and -Hhigh-imp- 

To motivate this step, consider the on-site Green func- 
tion of the Hubbard model G(w). Via the DMFT self- 
consistency condition, G(ui) defines an effective SIAM. 



For U t—> U c one indeed finds (e. g. within the NRG) that 
it makes no significant difference for the low-energy part 
of the solution of the resulting SIAM whether the full 
G(u>) is considered or the Green function with the Hub- 
bard bands removed. This means that in the iterative 
solution of the DMFT equations, the low-energy peak 
basically reproduces itself, and the high-energy degrees 
of freedom are rather unimportant. 

Alternatively, the first approximation can be charac- 
terized as follows: Let us (for a moment) look at the 
insulating solution for U — U c . Here the low-energy de- 
grees of freedom are absent (-ffiow-imp, ^low = 0) and the 
approximation reads: H high + Hhigh-imp + #im P >-> H ilap . 
The impurity spectral function for the left-hand side is 
given by two Hubbard bands centered at ±C//2 while it 
is given by two (5-peaks at ±E//2 for the right-hand side. 
So we can state that in step (i) of the approximation the 
finite bandwidth of the Hubbard peaks is neglected. 

(ii) For U i— ► U c , the width of the quasiparticle peak 
vanishes. This fact is used for the second approximation: 
We assume that in this limit it is sufficient to describe 
the low-energy degrees of freedom by a single conduction- 
band level, i. e. H\ OVJ i— > J2 a € cc\ r c tJ . Thereby we disre- 
gard any internal structure of the quasiparticle peak for 
U h- > U c . Equivalently, this means that the hybridization 
function can be represented by a single pole at u> — 0: 

ah = ■ ( 3 ) 

CO 

A given hybridization function A (a;) fixes the effective 
impurity problem. Due to 

A M = J2 7^ V - ( 4 ) 

V w - (e fc - /x) 

the one-pole structure of A(u>) corresponds to an n s = 2 
site single-impurity Anderson model (SIAM): 

#2-site = tdd\d a + Ud\d^d\di 

a 

+ 5> c t CCT + E y ( d K + hx -) (5) 

with the hybridization strength V = vSjv- 

Combining (i) and (ii) we obtain: -ffsiAM ^ -^imp + 
-ffiow + H iow-imp | — * ^2-site- With these two approxima- 
tions we can run through the DMFT self-consistency cy- 
cle. In the one-pole ansatz (^) for the hybridization func- 
tion, Ajv is the weight of the pole. The index TV refers 
to the AT-th. step in the iterative solution. Our goal is to 
calculate Ajv+i- 

We restrict ourselves to the manifest particle-hole sym- 
metric case. The chemical potential is set to // = U/2. 
The on-site energies in (^) are thus given by = tn — 
and e c = U /2. For the hybridization strength we have 
V = \/Ajv ' * as U i * U c . The two-site impurity 
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model is simple enough to be solved analytically |18|,|19|. 
For small V there are two peaks in the impurity spectral 
function atwsi ±U/2 as well as two peaks near to = 
which can be considered as corresponding to the quasi- 
particle resonance of the infinite (n s — oo) system. The 
weight of this "resonance" can be read off from the ex- 
act solution |ig| ]; up to second order in V/U and for the 
particle-hole symmetric case it is given by: 



l&V 2 



36 

U 2 



(6) 



In the self-consistent solution, z is also the quasiparti- 
cle weight which determines the low-energy behavior of 
the (local) self-energy of the lattice problem: 



S(w) = £7/2 + (1 - z- 1 )^ + 0(iv 2 ) 



(7) 



For a homogeneous lattice and a local self-energy the on- 
site Green function of the Hubbard model can be written 



as: 



G{uj) = / de 



uj — (e — /i) — T,(uj) ' 



(8) 



where p(e) is the free (U = 0) density of states. Using 
cq. ([?]) we obtain: 



G(uj) = z f de-^- + G {incoh -\uj) , 



u — ze 



(9) 



where the first part represents the coherent part of the 
Green function (G^ coh ''(uj)), and the second (incoherent) 
part can be disregarded for small excitation energies to <— > 
0. 

The integration can formally be carried out by means 
of a continued-fraction expansion which for a symmetric 
density of states p(e) reads: 

G^(oj) ee z [ de-^- = G^iz-'co) 
J uj — ze 

= l/(z- 1 cu-bl/(z- 1 u J ~b 2 2 / •■■)). (10) 



The expansion coefficients b n are related to the moments 
M n of the U = density of states. The first coefficient 
bi is given by: 



b\ = M 2 = J de e 2 p{e) 



(11) 



The second moment Mi is easily calculated by evaluating 
an (anti-)commutator of the form ([ [[c, H ]-, H ]_, c']+) 
which yields: 



M 2 = J2 i 



We thus obtain: 

G( coh -\u)) = 



uj - z 2 M 2 F(lu) ' 



(12) 



(13) 



where we have F(us) — \/ui + 0(ui 2 ) for the remainder. 

Starting from eq. (§) in the iV-th step, the DMFT self- 
consistency equation, 



A( W ) 



(e d — fi) — EH - G(uY 



(14) 



yields a new hybridization function A(uj) for the (iV+1)- 
th step. With eqs. (0) and Q we get: 



A(cj) = zM 2 F(uj) 



(15) 



for low frequencies u> t- 
structure, 



0. Insisting on the one-pole 



(16) 



for U i— » U c , we must have F (lu) = 1/lu. This amounts to 
replacing the coherent part of the on-site Green function 
by the simplest Green function with the same moments 
up to the second one. 

From eqs. (||) and (jlj) we thus have: 



A 



JV+l 



36 

U 2 



Mo A 



N 



(17) 



The coefficient of the (N + l)-th iteration step is thereby 
expressed in terms of the coefficient of the TV-th step. 
This is our main result. For U — U c the DMFT equa- 
tions are linearized, they are reduced to a simple linear 
algebraic equation which determines the evolution of a 
single parameter (A at) under subsequent iterations. 

The linearized mean-field equation ( |l7| ) has only one 
non-trivial solution with A^v+i = A^v which occurs for 
U = U C with U 2 = 36M 2 . Any U <U C gives A N+X /A N > 
1, so that An increases exponentially with iteration num- 
ber. This indicates the breakdown of the one-pole ap- 
proximation. For any U > U c the weight Ajv decreases 
exponentially with increasing iteration number. This cor- 
responds to the vanishing of the quasiparticle peak in the 
insulating regime. Consequently, U c has the meaning of 
the critical interaction for the Mott-Hubbard transition 
and its value is given by: 



U c = 6 



IJdEe*p(e)=6 j£t 



(18) 



For a lattice with nearest-neighbor coordination num- 
ber q and hopping integral t = \tij\ between nearest- 
neighbors i and j, we have: U c — 6t y/q. 

The result (|l|) has been derived within the DMFT 
which becomes exact in the limit q <— > 00. With the 
usual scaling for the hopping integral, t — t* / 1 ^fq and 
t* = const, km, we have: 



U c = 6f 



(19) 



However, eq. (18) may also be used for finite-dimensional 
systems where the DMFT is effectively the approxima- 
tion of a purely local self-energy functional. 
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III. DISCUSSION 
A. Bethe and hypercubic lattice 

For the Bethe lattice with infinite coordination num- 
ber and scaling t — t* / \/q, the free bandwidth is given by 
W = U*. So we expect from eq. @ the Mott-Hubbard 
metal-insulator transition to occur at U c — 1.5W. This 
result is in very good agreement with the result from 
the Projective Self-consistent Method (PSCM) [0JlO| 
£^c,pscm~ 1.46 W and with recent calculations using the 
Numerical Renormalization Group (NRG) method @] 
which yield J7 c ,nrg ~ 1-47W. It also agrees well with 
the value of U c ~ 1.511^ obtained in the NRG calcu- 
lations of Shimizu and Sakai |2(]]. The earlier IPT re- 
sult C/ c .ipt ~ 1.65W jl(J overestimates the critical U as 
compared to the other, non-perturbative methods. The 
Random Dispersion Approximation (RDA) |L3| predicts 
a considerably lower critical value [/ Cj rda ~ W. The ori- 
gin of this discrepancy, however, is presently not clear. 

On the infinite-dimensional hypercubic lattice with the 
scaling t = t* / y/q, we expect the metal- insulator transi- 
tion to occur at U c = 6t*. Again, this agrees well with 
the NRG calculations @ where the value U c « 5.80i* 
has been found. 

The existence of a metal-insulator transition in the hy- 
percubic lattice at a finite U c is not at all clear, consider- 
ing the fact that the free density of states is Gaussian, i.e. 
has no cutoff. In any case, the actual bandwidth (which 
is infinite for a Gaussian density of states) cannot play a 
role for the value of U c . It is much more plausible that 
it is the effective bandwidth (which is proportional to 

\J J dep(e)e 2 ) that has to be taken as a measure for U c . 

Our analysis also shows that U c is roughly independent 
of the details of the lattice structure and only depends 
on the local quantity J2j Hj ■ This result can quite natu- 
rally be understood when the electrons are considered as 
getting localized at the transition. In this case the elec- 
trons would only see their immediate surrounding which 
is the same for both the infinite dimensional Bethe and 
the hypercubic lattice. 

B. Two-sublattice model 

Let us now work out the predictions of the linearized 
DMFT for inhomogeneous lattices, i. e. lattices with re- 
duced (translational) symmetries. The presumably sim- 
plest but non-trivial case is a Bethe lattice that consists 
of two non-equivalent sublattices Qi and Q 2 where each 
site of Q a has q a nearest neighbors that belong to 
(with a — 2 for a = 1 and a = 1 for a = 2). We consider 
the limit of infinite coordination numbers q\ , qi i— ► oo 
with < qi/qs < oo. As for the homogeneous case 
qi = q2, it can be shown that the Hubbard model on 



the inhomogeneous lattice remains well-defined and non- 
trivial if the hopping integral is scaled appropriately, e. g. 
t = t*/^qi+q 2 = t**/y/qj with t*,t** = const. As a 
consequence, the self-energy T, a (ui) is local but sublattice 
dependent. The lattice problem can be mapped onto two 
impurity models that are characterized by hybridization 
functions A a (uj). The DMFT self-consistency equations 
read: 

A a (w) =u-[e d -ix)- E a (w) - Gaicj)- 1 , (20) 

where G a (tu) is the on-site Green function for a site i 
within the sublattice a. One easily verifies that the free 
(U = 0) local density of states on each sublattice is sym- 
metric and that /i = [7/2 at half-filling. Furthermore, 
with qi , q 2 i— > oo , we obtain from the lattice Dyson equa- 
tion: 

G^lu)- 1 = UJ + [1 - E a (w) - q a t 2 Gar(iu) . (21) 

The linearized DMFT for U i— ► U c iterates the one-pole 
ansatz A a (ui) = /ui. From eqs. ( p0[ ) and ( pl| ) we have 
A a (w) = q a t 2 Ga(oj) for a — 1,2. This implies that the 
quasiparticle peak for the sublattice a with weight z„ = 
(36/U 2 )A^ generates a corresponding peak in A a (uj) 
with the weight A^L = q a t 2 Za- We thus get: 

^N+i=J2 K ^ A N^ ( 22 ) 



where the 2x2 matrix K(Z7) is defined as: 

A fixed point of K(C7) corresponds to a self-consistent 
solution. Let X r (U) denote the eigenvalues of K(J7). 
We can distinguish between two cases: If |A r (Z7)l < 1 
for r = 1 and r = 2, there is the trivial solution 
limN^oo A^ = only (insulating solution for U > U c ). 

On the other hand, if there is at least one X r (U) > 1, A^ 
diverges exponentially as N i— > oo (metallic solution for 
U < U c ). The critical interaction is thus determined via 
the maximum eigenvalue by the condition: 

A ma x(C/ c ) - 1 • (24) 

This yields: 

U c - 6t tfq-^- 2 , (25) 

i. e. the geometrical mean of the critical interactions of 
two homogeneous Bethe lattices with coordination num- 
bers q\ and q 2 , respectively. The result recovers the ho- 
mogeneous case qi = q 2 and correctly gives U c = for 
the atomic limit q\ i— ► or q 2 i— ► 0. 

The analysis can be generalized straightforwardly to 
an arbitrary number of s sublattices Q%, . . . , Q s . We con- 
sider a Bethe lattice where each site of the sublattice Q a 
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has (q a — 1) nearest neighbors belonging to the sublat- 
tice Q a+ and one nearest neighbor in the sublattice Q a _ 
where a± = a ± 1 except for a = s (here a + = 1) and 
a = 1 (a_ = s). In the limit q a i— > oo with fixed pairwise 
ratios < q a /qp < oo, we have: 

Gq^w) -1 = w + ^ - - q a t 2 G a+ (u>) , (26) 

and the argument proceeds as above. We finally arrive 
at the mean-field equation ( p2j ) with K(£7) being an s- 
dimensional matrix with s non-zero elements: 



Let G^^w) = G a {uj) be the on-site Green function 



K(U) 



This implies: 



36 1 2 
U 2 



( o 



q 2 




9s- 1 
/ 



(27) 



l/2s 



£/c = 6i q a 



(28) 



\ct=l 



Again, this is plausible since q a = for any a would mean 
to cut the lattice into unconnected pieces of finite size, 
and the Mott transition becomes impossible (U c = 0). 

Also the s i— > oo limit of eq. (^8|) is meaningful: Con- 
sider e. g. q a =i 7^ q = <?2 = 93 = • • •• This describes 
a Bethe lattice with coordination number q for all sites 
except for one distinguished impurity site with coordina- 
tion number q±. As expected physically, U c is unaffected 
by the presence of the impurity. Furthermore, in any 
case where one changes the number of nearest neighbors 
of a finite number of sites only, the value for U c remains 
unchanged. 



C. General inhomogeneous Bethe lattice 

We finally tackle the "inverse" problem: Given a ma- 
trix K, is there a realization of a (Bethe) lattice such 
that the critical interaction is determined by the max- 
imum eigenvalue of K? For this purpose we consider 
the Hubbard model with nearest-neighbor hopping on a 
general inhomogeneous Bethe lattice where each site i 
may have a different coordination number. Remaining 
spatial symmetries are accounted for by classifying the 
lattice sites into sublattices Q a that consist of equiva- 
lent sites only. By q a p we denote the number of nearest 
neighbors of a site i £ Q a that belong to the sublattice 
Qp. We are interested in the limit q a p <— > oo with < 
qup/q-yS < oo since this implies a local but a-dependent 
self-energy S a (w). Within the DMFT this Hubbard 
model is mapped onto impurity models which are labeled 
by the sublattice index a. The self-consistency conditions 
are given by eq. (pc| ) . 



for a site i in Q a , and G, 



(n) 



••aia(w) = ((c iCT ;ct CT )) w 

the off-site Green function for n-th nearest-neighbor sites 
* £ Qa n and j £ Q a . G^ depends on the sublattice 
indices that are met along the (unique) path from j to i. 
Via its equation of motion, Gq ^ couples to the nearest- 
neighbor off-site Green function Ga}a- 

(w + M - S„(«))GW (w) = 1 + tJ2 Q aai G£U u ) ■ ( 29 ) 



For a Bethe lattice the nearest-neighbor Green function 
Gala can only couple to Ga}a ia and to G^ again. More 
generally, the equation of motion for the n-th nearest- 
neighbor off-site Green function reads: 



(U> + fl - Y, an {uj))G a a }...a l a{^) = tG K a n _\ '... aia {uj) 

+t q an a n + 1 G a J^ 1 \. aia (uj) , (30) 

where in the second term on the r.h.s. we have used the 
approximation q a p -1» q a p which becomes exact in the 
limit of infinite coordination numbers. The infinite series 
defined by eqs. (|9|) and ( |30| ) can formally be summed 
up. This yields: 



G a (uj) 1 = LP + fl - ^ a (u) - ^y^^qapGpju) 

P 



(31) 



The mean-field equation of the linearized DMFT thus 



has again the form (22) where the if-matrix is given by: 



tr (tt\ m2 
K a p(U) = -jj^-q a p ■ 



(32) 



The critical interaction is given by 6i times the maximum 
eigenvalue of the coordination-number matrix q. Note 
that for a general (non-symmetric), irreducible matrix 
with non-negative elements, the eigenvalue with maxi- 
mum absolute value is real and non-negative (Perron- 
Frobenius theorem We thus conclude that any 

quadratic matrix K with non-negative elements can be 
related to the Mott transition on a Bethe lattice with 
certain (infinite) coordination numbers. 



D. Hypercubic films 

As a more realistic example for the Mott transition on 
an inhomogeneous lattice we consider a hypercubic Hub- 
bard film. A D-dimensional film is built up from a num- 
ber d of (D — l)-dimensional "layers". For a hypercubic 
film these layers are cut out of the usual D-dimensional 
hypercubic lattice. A set of Miller indices [x%, x%, . . . , Xjj\ 
characterizes the film-surface normal direction. The most 
simple films are those with low-index surfaces given by 
x\ = ■ ■ ■ = x r — \ and x r +\ = ■ ■ ■ = xd = 0. For any site 
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in the film except for sites at the film surfaces there are 
qi = 2D — 2r nearest neighbors within the same layer and 
q 2 = r nearest neighbors in each of the adjacent layers; 
the total coordination number is q = q% + 2q 2 = 2D. 

For D i— ► oo the Hubbard model is well defined with 
the usual scaling of the hopping t = t*/y/2D, the self- 
energy becomes local but layer dependent and dynamical 
mean- field theory is exact [Q. The lattice problem is 
mapped onto a set of d impurity problems. The DMFT 
self-consistency conditions are given by eq. (^) where 
the index a now has to be interpreted as the layer index: 
a = 1, . . . , d. 

The linearized DMFT can be developed as in Sec. II. 
Equation (||), however, is no longer valid and must be re- 
placed by the Dyson equation corresponding to the given 
film geometry. The coherent part of the on-site Green 
function for a site in the layer a is given by: 



G a (u) = z a G a (u) = 



u-M^ a) F a (u) 



(33) 



where z a is the layer-dependent quasiparticle weight, 
z a = (1 — dT, a (iO + )/dui)~ 1 and G a (uj) is the on-site 
element of the free (U — 0) Green function but calcu- 



lated for the renormalized hopping t{ 



with Zi = z a for a site i in the layer a. In the ex- 
pression on the right, denotes the corresponding 



second moment which is calculated as M, 



z a (qiz 



Q2Z a -i+q2Z a +i)t 2 ■ This means that the linearized mean- 
field equation has again the form ( p2| ) with the following 
e?-dimensional tridiagonal matrix: 



K(U) 



36 t 2 
U 2 



\ 



qi <i2 
<?2 qi <?2 

<?2 qi 



J 



(34) 



Its eigenvalues are the zeros of the d-th degree Chebyshev 
polynomial of the second kind p3fl . From the maximum 
eigenvalue we obtain: 



U c = 6t \ qi + 2q 2 cos 



1 



(35) 



Equation (|35|) describes the complete thickness and ge- 
ometry dependence of the critical interaction for the 
Mott-Hubbard transition in hypercubic Hubbard films. 
In the limit of thick films d <— > oo one recovers the bulk 



value U c = 6t^qi + 2q 2 . For d < oo the critical inter- 
action not only depends on the film thickness d but also 
on the geometry of the film surface which is character- 
ized by r. Varying r we can pass continuously from the 
most closed (r — 1) to the most open (r = D) surface 
geometry. For r = 1, i. e. a (1000...) film surface, a site 
in the topmost layer has qs = q\ + q 2 = 2D — 1 near- 
est neighbors to be compared with q = 2D in the bulk. 
For D i * oo the local environment of the surface sites 



is essentially the same as in the bulk, i. e. surface effects 
become meaningless. Consequently, we get U c = 6t* , i. e. 
the bulk value irrespective of d. For r = D one obtains 
the open (1111...) film surface. The surface coordination 
number is reduced to q§ — q 2 — D. The critical inter- 
action is U c = 6t* cos(ir / (d + 1)) which is smaller than 
6t* for any d. 

Equation ( ]35| ) can also be applied to finite-dimensional 
films (D < oo) if one additionally assumes the local ap- 
proximation for the self-energy functional to hold. For 
D = 3 simple-cubic films with a thickness ranging from 
d = 1 up to d = 8 and for sc(100), sc(110) and sc(lll) 
film surfaces, the prediction ( |3q ) of the linearized DMFT 
has been tested in Rcf. \22\ by comparing with the re- 
sults for U c of a fully numerical evaluation of the DMFT 
equations using the exact diagonalization of small im- 
purity models (n s =8). It is found that the linearized 
DMFT qualitatively and - as far as can be judged from 
the numerical evaluation - also quantitatively predicts 
the correct geometry and thickness dependence of U c [E2| . 



E. Critical exponent and critical profiles 

So far the discussion was restricted to the calculation 
of the critical value U c which is derived from a linear 
homogeneous mean- field equation [eqs. (|l7| ) and (p2|)j. 
To determine the critical behaviour of the quasiparticle 
weight z for U i— > U c , one has to go beyond the linearized 
DMFT. For this purpose a simple generalization of the 
arguments in Sec. II is necessary. 

We replace the second-order result for the quasiparticle 
weight z = 36V 2 /U 2 [eq. (0)] by the result up to fourth 
order in V/U: 

V 2 ( V 2 \ 

- = 36 ^( v 1 - 4 Vj • (36) 

With the same steps as before, one arrives at: 

1 - J j^M 2 A N , (37) 

which is a non- linear equation for the "mean field" A. 
The self-consistency requires A^+i = = A. Solving 
for A yields: 



A = itfc (U c - U) 



(38) 



where we have already expanded the right hand side in 
powers of (U c —U), The result for the quasiparticle weight 
near U r is 



18 U C -U 
TT U c 



(39) 



This equation is, of course, only valid for U < U c . We 
obtain a linear vanishing of the quasiparticle weight near 



G 



the metal-insulator transition as in the Brinkman-Rice 
approach @ and the PSCM @. 

For the Hubbard model on an inhomogeneous lattice, 
the self-energy and thus the quasiparticle weight is site 
or sublattice dependent: z a = (1 — dY, a (iO + )/du})~ 1 . 
Within the linearized DMFT [eq. @], the critical 
interaction is determined from the largest eigenvalue 
A max (J7 c ) = 1 in the eigenvalue problem 



Z0 



(40) 



The corresponding eigenvector z a — z a (U c ) describes the 
critical profile of the quasiparticle weight. The profile 
is uniquely determined up to a normalization constant. 
While z a (U) h- > for each a as U i— > U c , the ratios 
z a {U c )/ zp(U c ) remain to be non-trivial. For example, in 
the two-sublattice model characterized by e q. (|23| ), the 
critical profile is given by zi(U c ) / Z2{U C ) = \fq\Jq-2- 

To determine the critical behaviour of the a-dependent 
quasiparticle weight for U 1— > U c but U < U c , we again 
have to expand up to fourth order in V a /U. This yields 
the following mean-field equation: 



z a (U)=Y,K a p{U) z p {U) 
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4(u) 



For the two-sublattice model we obtain: 
36 JqH. U c -U 



For qi — q 2 this result reduces to eq. (| 



(41) 



(42) 



IV. SUMMARY 



transition, such as the existence or absence of a hystere- 
sis, the order of the transition, etc. Its advantage is that 
it can be easily generalized to a variety of geometries. 
To obtain the critical interaction, it is sufficient to find 
the maximum eigenvalue of the respective coordination- 
number matrix, the dimension of which is determined 
by the remaining spatial symmetries (see e.g. eq. (|32"|)). 
The analytical results for e. g. the metal-insulator tran- 
sition in thin Hubbard films have been checked against 
numerical solutions of the full DMFT equations, and the 
geometry dependence has been found to be essentially 
the same, in the DMFT and in the linearized DMFT. 

The main approximation which our approach is based 
on, is the one-pole structure of the effective hybridization 
function. It is obvious that the linearized DMFT can in 
principle be extended by taking into account more states 
for the effective conduction band (in the quasiparticle 
peak and/or the Hubbard bands). However, its main 
advantage - the possibility to obtain analytical results 
for U c - would then be immediately lost. 

It would be very interesting to sec whether experi- 
ments on Mott-Hubbard systems in different geometries 
will show similar trends as predicted by the linearized 
DMFT. 
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We have discussed a "linearized" version of the Dy- 
namical Mean Field Theory which allows for the analyt- 
ical calculation of the critical interaction for the Mott- 
Hubbard metal- insulator transition at T = 0. The main 
result is: 



U c = 6^jjd66^p(s)=6 ] JTA' ( 43 ) 

which shows that it is the second moment of the non- 
interacting density of states which determines U c . This 
is contrary to the Hubbard-III approximation where the 
critical interaction is determined by the free bandwidth 
W and to the Brinkman-Rice approach where U c is given 
in terms of the U — kinetic energy eg. The values for U c 
obtained with the linearized DMFT have been compared 
with the available results from numerical solutions of the 
full DMFT equations, and a good agreement is found. 

The linearized DMFT is of course not able to answer 
detailed questions about the nature of the Mott-Hubbard 
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